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Neutrino mixing anarchy is the hypothesis that the leptonic mixing matrix can be described as 
the result of a random draw from an unbiased distribution of unitary three-by-three matrices. In 
light of the recent very strong evidence for a nonzero sin 2 2(9i3, we show that the anarchy hypothesis 
is consistent with the choice made by the Nature - the probability of a more unusual choice is 
44%. We revisit anarchy's ability to make predictions, concentrating on correlations - or lack 
thereof - among the different neutrino mixing parameters, especially sin 2 #13 and sin 2 #23- We also 
comment on anarchical expectations regarding the magnitude of CP-violation in the lepton sector, 
and potential connections to underlying flavor models or the landscape. 
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The flavor puzzle has befuddled generations of particle physicists. Since the first years of the quark model and the 
first successful description of flavor-violating weak processes, the pattern of fermion masses and mixing parameters 
seems to hint at the existence of some yet-to-be-uncovered organizing principle. The main idea is that new hidden 
I symmetries - global or local, spontaneously or explicitly broken - "explain" the fact that the charged-fermion masses 
are very hierarchical and that the quark mixing matrix is very close to the identity matrix. 

The discovery of nonzero neutrino masses and lepton mixing in the end of the last century added new pieces to the 
flavor puzzle. In particular, the structure of U, the leptonic mixing matrix [22| . seems to be providing qualitatively 
different information. Unlike the quark mixing matrix, U cannot be understood as an identity matrix "perturbed" 
Cl^ by small, hierarchical, off-diagonal matrices. Qualitatively speaking, all elements of the leptonic mixing matrix are 
large: \U ai \ = 0(l). 

The flavor literature is densely populated with ingenious attempts to identify the organizing principle behind U. 
The simplest idea, arguably, is to postulate that there is, in fact, no organizing principle behind U. Neutrino mixing 
anarchy [l| is the hypothesis that U can be described as the result of a random draw from an unbiased distribution 
of unitary 3x3 matrices, as discussed in detail in Q. Several years ago 0, we proposed that a Kolmogorov-Smirnov 
(KS) test could be employed in order to test the anarchy hypothesis. At the time, anarchy provided a very good fit to 
[ the neutrino oscillation data and we further proposed to use the KS statistic in order to make predictions concerning 
yet-to-be-determined mixing parameters assuming that the anarchy hypothesis is correct. In particular we predicted 
that, at the two sigma level, |£/ e 3| 2 > 0.011. Here we revisit the anarchy hypothesis in light of new experimental 
developments regarding the leptonic mixing matrix. 
We parameterize U as, 
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where A3 = diag(l,— 1,0) and As = diag(l, 1, — 2)/V3 are Gell-Mann matrices, and — sin&y, cy = cosOij, 
ij — 12, 13, 23. (pi 2 are unphysical phases that can be absorbed by rephasing the charged lepton fields, Xi,2 are 
potentially physical "Majorana" phases, and S is the "Dirac" CP-odd phase that is potentially observable in long- 
baseline neutrino oscillation experiments. 

In 2011, data from T2K [|, MINOS @, and Double Chooz Q were consistent with \U e3 \ 2 = sin 2 6 13 ~ 0.02, 
and ruled out, for the first time, \U e3 \ 2 — at the three sigma level. These data were in agreement with previous 
hints, from combinations of different neutrino experiments, that |J7 e3 | 2 did not vanish Q. Very recently, Daya Bay 
released the analysis of their first data and claimed five sigma evidence that \U e 3\ 2 is not zero [!|. Its data point to 
\U e 3\ 2 = 0.023 ± 0.004, quite consistent with all previous hints for a nonzero |£/ e 3| 2 - A few weeks after [§], RENO 
released its first results [10| . which are in agreement with those from Daya Bay and rule out | t/ e 3 1 2 = at more than 
six sigma. The RENO results point to | C/ e 3 1 2 = 0.026 ± 0.004. The other two mixing angles that parameterize U are 
relatively well-measured from the world's neutrino data (for recent global analyses, see [ill HI])) an d both are large. 

The current data can be summarized as [H, [l2[ 

sin 2 6>i2 = 0.312±°;^, sin 2 6> 23 = 0.52 ± 0.06, sin 2 13 = 0.023 ± 0.004, (2) 

where the last measurement is from Daya Bay. Henceforth, instead of attempting to combine all measuements of 
sin 2 0i3, we will use the Daya Bay result when it comes to the experimental value of sin 2 #i 3 . This is both a very 
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good approximation to the world neutrino data and, for our purposes, turns out to be a more conservative choice 
(slightly smaller central value and larger uncertainty) . 5 and the potentially physical Majorana phases remain virtually 
unconstrained. 

The anarchy hypothesis requires the probability measure of the neutrino mixing matrix to be invariant under 
changes of basis for the three generations. It immediately leads to a unique probability distribution for the mixing 
angles. These can be read off from the invariant integration Haar measure Q, 

dU = ds\ 2 A dcf 3 A ds 2 3 A dS A dr) A d<j)\ A dfo A dxi A d\2, (3) 

up to an overall normalization factor. It yields, for example, probability distributions in sin 2 20 that are the same for 
012, Gig, and 23 0. 

In [3j, we derived the probability that the anarchy hypothesis is consistent with experimental data as 

P 3 KS = £3 (l - log e 3 + \ log 2 e 3 ) , (4) 

where 

e 3 = 2 sin 2 8 12 x 2 min(sin 2 23 , cos 2 23 ) x 2 (l - cos 4 0i 3 ) . (5) 

The expression reflects the fact that the Haar measure Eq. @ is flat in s 2 2 , s| 3 , and c 4 3 . Eq. Q is a function of the 
three mixing angles and is computed by marginalizing over all other parameters, the phases <5, Xi,2, 4*1,2 in Eq. ([TJ. In 
Eq. (O, 'min' indicates that, depending on whether 23 is in the "light side" (02 3 < 7r/4) or the "dark side" (023 > 7r/4) 
[H of its parameter space, one should use either sin 2 23 or cos 2 2 3, respectively. 

For the central values of the mixing parameters in Eq. ([2]), e 3 = 0.054 and P^ s = 44%. P^ s can be interpreted as 
the probability that a random draw for the mixing matrix will yield a result that is more "unusual" than the one made 
by Nature. In a nutshell, the anarchy hypothesis is consistent with our current understanding of lepton mixing, and 
the observed values for the mixing parameters are quite typical of the ones obtained by randomly drawing a mixing 
matrix from an unbiased distribution of unitary 3x3 matrices. 

Fig. H] depicts P^ s as a function of sin 2 0i 3 , for values of sin 2 012,23 that independently span their respective allowed 
regions at the three-sigma confidence level [l2| ■ The dotted vertical lines indicate the values of sin 2 0i3 allowed by 
Daya Bay at the three-sigma confidence level and allow one to conclude that, within the three-sigma experimentally 
allowed region of the parameter space, P^ s > 24%. 




FIG. 1: P^ s as a function of sin 2 6\z, for values of sin 2 612,23 that independently span the respective allowed regions at the 
three-sigma confidence level [HI]. The solid [blue] line corresponds to the best-fit values, while the dashed [red] lines define the 
largest and smallest allowed values of P^ s ■ The [green] dotted vertical lines indicate the values of sin 2 #13 allowed by Daya 
Bay at the three-sigma confidence level. The best estimate for (sin 2 #13, P 3 XS ) = (0.023,44%) is indicated by the black dot. 
The horizontal line segments indicate different experimental measurements of sin 2 #13 (one-sigma error bars) . In the case of 
MINOS and T2K, we assume a normal neutrino mass hierarchy. 



In order to see how the measurements of the three mixing angles contribute to such an excellent agreement with the 
anarchy hypothesis, we can define a KS probability P 2 S &jk), hj, fc = 1,2,3 as a function of two mixing angles 
by marginalizing over the "other" mixing angle. Following the procedure discussed in detail in [3], 



^ S (%, J fe) = £2(l-loge 2 ), 



(6) 
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where 



£2(6*12, 6*23) = 2 sin 2 6»i2 x 2 min(sin 2 23 , cos 2 2 3 
£2(^12, #13) 
£2(6*13, 623) 



2sin 2 6»i 2 x 2(l-cos 4 6»i 3 ), 
2 min(sin 2 9 23 , cos 2 #23) x 2 (l - cos 4 0ig 



(7) 
(8) 
(9) 



Fig. [2] depicts contours of constant P 2 KS {9ijsk) = 68.3% ("one sigma") and P2 KS (0ij,ik) = 95.5% ("two sigma") in 
the sin 2 0y x sin 2 9jk plane, for the different combinations of Lj, k. The plots also depict the best fits to the current 
neutrino data at the one and three sigma confidence levels [H, [Hj . Note that we assume that the data responsible for 
the allowed values of sin 2 9i 2 , sin 2 #13, and sin 2 6*23 are completely uncorrelated. While this is not completely correct, 
it is a good enough approximation for the intentions of this manuscript. The anarchical predictions for the pair- wise 
values of sin 2 9 are in good agreement with experimental data (the one-sigma "tension" in the sin 2 9i 2 x sin 2 #13 plane 
is clearly not statistically significant). This adds credence to the already advertised fact that the anarchy hypothesis 
is in very good agreement with the neutrino oscillation data. 
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FIG. 2: Constant P 2 (%ifc) = 68.3% ("one sigma," [black] dashed contours) and Pf {Oij,ik) = 95.5% ("two sigma" [red] 
solid contours) contours, in the three distinct sin 2 @ij x sin 2 Ojk planes. The shaded regions indicate the currently allowed region 
of the parameter space at the three sigma and one sigma (black ellipses) confidence levels, according to [3, Oil- Correlations 
among the experimentally allowed values of different sin 2 Oij were not taken into consideration. 



Our results do not, of course, imply that anarchy is necessarily correct. Strictly speaking, they simply mean that 
the neutrino data do not falsify the anarchy hypothesis. Flavor models also make predictions for the values of the 
mixing angles. We will not attempt to summarize all different possibilities identified in the literature, but will describe 
in some detail one concrete example, for illustrative purposes. 

Attempts to identify order in the leptonic mixing matrix follow different paths. One is to postulate that a new 
organizing principle, at leading order, leads to (i) |?7 T j| = \U^i\, "maximal" atmospheric mixing, and (ii) |£7 e 3| = 0. 
These combined translate into sin 2 6*23 = 1/2, sin #13 = 0. Other zeroth order predictions may also apply, including 
I C/ e 2 1 2 = |^2| 2 = \U t2 \ 2 = 1/3, mi = m 2 = 0, etc. Higher order effects will lead to deviations from these zeroth 
order predictions. These effects - their magnitudes and associated flavor structure - are, of course, model dependent. 
If, however, the higher order corrections are generic, e.g., flavor blind and governed by a single small parameter, 
one generically predicts that deviations of atmospheric mixing from maximal (which can be parameterized by the 
deviation of cos 26*23 from zero) are correlated with deviations of sin 2 #13 from zero. This is discussed in some detail 
in, e.g., [l3|. More quantitatively, one predicts, under the circumstances outlined above, 

sin 2 6»i 3 = Ccos 2 26>23 = C(l - 2 sin 2 9 23 ) 2 , (10) 

where C — 0(1), a proportionality constant that is either a free parameter or is specified by the model in question (for 
a concrete, recent example, see [lj]). This relation in the sin 2 #23 x sin 2 6*13 plane is depicted in Fig.|3]for C € [0.8, 1.2] 



4 



(parabolic [blue] region). We refer to this region of the parameter space as the prediction of the 'ordered hypothesis.' 
The figure also depicts the experimentally allowed values of sin 2 9 23, sin 2 9 13 at the one and three sigma levels, and 
the region of the parameter space preferred by anarchy at the one and two sigma levels, as in Fig. [2ftop-right). 
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FIG. 3: Prediction of the ordered hypothesis in the sin 2 823 X sin 2 #13 plane ([blue] parabolic contours), dictated by Eq. (|10|) 
for C £ [0.8, 1.2]. The light [yellow] curve corresponds to C = 1. The currently allowed region of the parameter space and the 
expectations from the anarchy hypothesis, both in Fig. [2jtop-right), are also depicted. 

Fig. [3] reveals that the ordered hypothesis and the anarchy one prefer somewhat different regions of the currently 
allowed sin 2 #23 x sin 2 #13 parameter space. The anarchy hypothesis does not strongly prefer any region of the 
experimentally available space. It does, however, favor maximal sin 2 #23 = 1/2 and "large" values of sin 2 # 13 . On 
the other hand, the ordered hypothesis, in light of the Day Bay result, rules out sin 2 #23 = 1/2, instead preferring 
cos 2#23 ~ ±0.1 (this point was recently also emphasized in [Hj]). It is also curious to note that C < 0.5 values are 
disfavored. 

Precision measurements of the neutrino oscillation parameters may ultimately favor anarchy versus order, or vice- 
versa. The values of the parameters are such that an improved determination of sin 2 #23 will provide the most 
discriminating power. If one interprets the width of the blue region in Fig. [3] as indicative of the uncertainty in the 
ordered predictions, next-generation experiments sensitive to i5(sin 2 #23) ~ 0.02 - an uncertainty of a few percent - 
would be required to qualitatively change our understanding of structure in lepton mixing. The NOf A experiment, 
for example, is aiming at measuring, from z/ M disappearance, sin 2 26*23 at the 0.4% level for sin 2 26*23 = 1 [3; which 
translates into an uncertainty of 0.03 for sin 2 #23 = 0.5. Similar, albeit slightly worse, precision is expected from T2K 
[l7j . The fact that #13 is large implies that — > v e searches at T2K and NOf A, combined with reactor measurements 
of v e disappearance, will allow one to directly measure sin 2 #23- The precision with which sin 2 #23 can be measured 
will be dominated by the precision with which T2K and NOz/A can measure sin 2 #23 sin 2 2fli3.[23j] which is expected 
to be markedly worse than the one advertised for sin 2 #23 from disappearance, above. Interesting information is 
also expected fro m p recision measurements of the atmospheric neutrinos at, for example, the INO experiment (see, 
for example, [l8l. 1 1 9j] . and references therein). 

Similar arguments can be made in the sin 2 #12 x sin 2 #13 and sin 2 Qyi x sin 2 #23 planes. The circumstances here, 
however, are different, sin 2 Qyi is already known at the few percent level. This implies that constraints on successful 
ordered scenarios are either very stringent and the associated "predictions" are very tight (e.g., sin 2 812 may almost 
uniquely determine the value of sin 2 #13 and sin 2 #23) or correlations are either absent or very weak. In the sin 2 #12 x 
sin 2 623 plane, the anarchical prediction works almost "too well," as the currently three-sigma experimentally allowed 
region is entirely contained deep in the one-sigma anarchy hypothesis prediction. It is quite unlikely that an ordered 
hypothesis will lead to a significantly better, statistically speaking, a posteriori agreement with the data. 

The next obvious target for neutrino oscillation experiments is the discovery of leptonic CP-invariance violation, 
whose magnitude is governed by the Dirac phase S. For example, for neutrinos propagating in vacuum, P(y^ — > 
— P(Pfj, —> v e ) oc sin 5. Since the Haar measure Eq. §5§ is fiat in 5, the probability distribution of sin 5 is peaked 
at sin 8 = ±1 the anarchy hypothesis implies that "large" leptonic CP-invariance violation is quite probable. 

If the neutrinos are Majorana fermions, the Majorana phases xi,2 in Eq. (TTJ) are physical observables. Similar 
to that of <5, their probability distributions are flat in Xl,2, respectively. Majorana phases are known to affect the 
magnitude of the neutrino exchange contribution to neutrinoless double-beta decay (Of/3/3), and it is interesting to ask 
whether the anarchy hypothesis has any impact on the expected rates for these rare nuclear processes. The answer, 
unfortunately, depends on the value of the lightest neutrino mass, which is both experimentally unknown and not 
addressed by the anarchy hypothesis, which concerns only mixing parameters. Nonetheless, we would like to advertise 
that, if the anarchy hypothesis is correct and neutrinos are Majorana fermions, it is quite unlikely that the rate for 
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Ov(3(3 decay is vanishingly small. 

For light neutrinos, the amplitude for Of/3/3 is proportional to 

mee = U ei™i = m 1 (e i ^ +i ^/^c 12 c 13 ) 2 + m 2 (e-^+^^c 13 s 12 ) 2 + m 3 {e'^ 2 ^ s 13 ) 2 , (11) 

i 

using Eq. ([T]). It is well-known that m ee may completely vanish for a normal neutrino mass hierarchy (mi < mi < m 3 ) 
for mi w 0.005 eV. Fig. @] depicts a histogram of \m ee \ values for mi = 0.005 eV obtained by varying \i,2 and S 
according to their anarchical probability distributions. The [red] curve indicates the probability that a given value of 
\m ee \ or larger is obtained. Tiny |m ee | values are allowed, but in the majority of the cases (94.5%) |m ee | > 0.001 eV. 
Even though this level is beyond the reach of current experimental efforts looking for 0^/3/3, it may be probed by 
future multiton-scale experiments (see, e.g., [2(| and references therein). 
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FIG. 4: Distribution of 160,000 |ra ee values (see Eq. (JTTJ) obtained by scanning over the values of the CP-odd phases 5, \i,2 
assuming their probability distributions are as prescribed by the anarchy hypothesis. We assume a normal neutrino mass 
hierarchy and m\ — 0.005 eV, while the mixing angles are fixed at their best fit values, Eq. J2J, along with the mass-squared 
differences [12] • The [red] curve indicates the probability that a given value of |m ee | or larger is obtained. 

The anarchy hypothesis, on the other hand, implies no correlation between the Dirac or Majorana phases with the 
resulting baryon asymmetry in standard leptogenesis [2l| . This implies that, if the anarchy hypothesis is correct, the 
direct experimental verification of leptogenesis is impossible. Nonetheless, the positive observation of CP-invariance 
violation in the lepton sector would still serve as a nontrivial plausibility test. 

We conclude by speculating on potential underlying theories beneath the anarchy hypothesis. Some belong to the 
category of flavor-symmetry models based on simple Z7(l)'s (as discussed, e.g., in [2j) where the three generations 
of lepton doublets are assigned identical flavor quantum numbers. In this case, the random numbers merely reflect 
our lack of understanding of the detailed predictions of the models, and the KS test discussed here offers a positive 
verification of the hypothesis of identical quantum numbers for all the lepton doublets. 

Another class of theories is related to the concept of a string-theory landscape, where our Universe is hypothesized 
to be only one among a vast number of possible Universes, each with different physics laws and parameters. Virtually 
all discussions of the landscape, especially when it comes to predictions - or lack thereof - rely on the "anthropic 
principle," where our own existence, or more generically that of galaxies, stars, etc, plays a role in determining which 
Universes observers are allowed to occupy and, hence, observe. "Environmental pressure" is then defined such that 
physical parameters, while predicted by the ultimate fundamental physics to be random with complicated probability 
distributions, are forced towards very special ranges of values and the Universe we inhabit, along with its physics laws 
and the values of the parameters, tends to be "on the edge" for intelligent life, or some appropriate facsimile, to exist. 
The bottom line is that the observable physics parameters take on what often appears to be ridiculously unlikely 
values in order to allow observers to exist, but the values are expected to be as mediocre as possible. If the values of 
the lepton mixing parameters are also a consequence of the landscape, the KS test presented here can be viewed as a 
direct, faithful test of the landscape hypothesis assuming - and the test seems to corroborate this! - leptonic mixing 
parameters have nothing to do with our existence (2lj . 
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